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Abstract. High moments of multiplicity distributions of conserved quantities are 
predicted to be sensitive to critical fluctuations. To understand the effect of the 
^ . non-critical physics backgrounds on the proposed observable, we have studied various 

moments of net-proton distributions with AMPT, Hijing, Therminator and UrQMD 
Qjl ' models, in which no QCD critical point physics is implemented. It is found that the 

5 , centrality evolution of various moments of net-proton distributions can be uniformly 

described by a superposition of emission sources. In addition, in the absence of critical 
phenomena, some moment products of net-proton distributions, related to the baryon 
I number susceptibilities in Lattice QCD calculations, are predicted to be constant as a 

' function of the collision centrality. We argue that a non-monotonic dependence of the 

, moment products as a function of the beam energy may be used to locate the QCD 

' critical point 

(N ■ 

O 

o . 

1. Introduction 

I Heavy-ion reactions at high energy allow us to study the QCD phase diagram 
^ ' experimentally [Ij. At vanishing baryon chemical potential [fiB = 0), Lattice QCD 
calculations predict that a cross-over from the hadronic phase to the Quark Gluon 
Plasma (QGP) phase will occur above a critical temperature. The temperature range 
for the cross-over has been estimated to be about 170 - 190 MeV [2]. QCD based model 
calculations indicate that at large fiB the transition from the hadronic phase to the 
QGP phase could be of first order with a critical point at the boundary to the cross- 
over, the QCD Critical Point (QCP) [3]. The location of the QCP or even its existence 
are not confirmed {3]. The possibility of the existence of the QCP has motivated our 
interest to search for it with the RHIC beam energy scan program [5]. By decreasing 
the collision energy down to a center of mass energy of 5 GeV we will be able to vary 
the baryo-chemical potential from ~ to of about 500 MeV. 

A characteristic feature of a critical point is the increase and divergence of the 
correlation length ( ^ ) and of critical fluctuations. In heavy-ion reactions, finite size 
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effects, rapid expansion, and critical slowing down could wash out those effects. For 
example, the critical correlation length in heavy-ion collisions is expected to be about 2- 
3 fm [6] . A clear signature of a critical point in an energy scan would be non-monotonic 
behavior of non-gaussian multiplicity fluctuations. 

Recently, theoretical calculations have shown that high moments of multiplicity 
distributions of conserved quantities, such as net-baryon, net-charge, and net- 
strangeness, are sensitive to the correlation length ^ [7]. 

In Lattice QCD calculation with /i^ = 0, higher order susceptibilities of the 
baryon number, which can be related to the higher order moments of the net-baryon 
multiplicity distributions, show a non-monotonic behavior near Tc |8| . A similar behavior 
is expected for the finite /i^ region. Experimentally, it is hard to measure the net- 
baryon number while the net-proton number is measurable. Theoretical calculations 
show that fluctuations of the net-proton number can be used to infer the net-baryon 
number fluctuations at the QCP [9]. 

In this paper we study the energy dependence of net-proton multiplicity 
distributions for several models in terms of Skewness ( S ) and Kurtosis { k ). This 
is a feasibility study for the future data analysis from the energy scan at RHIC. 



2. Observables 



We introduce various moment definitions of the event-by-event multiplicity distributions: 
Mean, M = < N >, Variance, = < (AA^)^ >, Skewness, S = < {ANf >/a^ and 
Kurtosis, K = < (AN)* >/cr^ — 3, where AA^ = A^— < N >. Skewness and Kurtosis 
are used to characterize the asymmetry and peakness of the multiplicity distributions, 
respectively. They are also used to demonstrate the non-Gaussian fluctuation feature 
near the QCP, in particular a sign change of the skewness may be a hint of crossing the 
phase boundary [71 |T0] . 

To understand the centrality evolution of these moments, we introduce the Identical 
Independent Emission Source (IIES) assumption. Here the colliding system consists of a 
large number of emission sources and the final multiplicity of particles is the sum of the 
multiplicities from individual emission sources. The relation between various moments 
and the number of emission sources for the i^^ centrality can be expressed as: 

n n 

M- a' 
EN, Ni 

i=l i=l 
n n 

N- AT. 
(3) : TTT^ = = S\x), (4) : — ^ = = k{x) 

■t=l i=l 

where Mj, a,, Si,Ki {i = 1,2, ...n) are the moments extracted from the multiplicity 
distribution of the i^^ centrality and Aj is the corresponding number of emission sources. 
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M{x),a{x), S{x), K,{x) are various moments of the multiplicity distributions for each 
emission source. From Equs .(l)-(4), we obtain: 

1/Sf 



Mi 

(5) : — ^ 



i=l i=l i=l i=l i=l 

which shows the connection between the emission source distributions and the various 
moments of multiplicity distributions. To investigate the centrality evolution of those 

n 

moments, we fit the normalized mean value Mi/E with a function /(< Npart >)? 

i=l 

where < Npart > is the average number of participants. Then, we obtain: 



(6) : M(< Npart >] 

(7) : a{< Npart >) 

(8) : S{< Npart >) 



Y.MA*f{< Npart >] 



i=l 



\ i=i 



^Y^l/Sf)* fi< Npart >] 
\ i=l 



It 

(9) : «:(< Npart >) = 1/[{Y.I/k,) * /(< Npart >)] 



i=l 



Consequently, the centrality evolution of various moments can be uniformly 
described by the function /(< Npart >)• From those equations it also follows that 
the moment products, Sa, Ka/S and kct^, are constant as a function of < Npart >■ 



3. Results 



We calculated the various moments of net-proton {Ap = Np — Np) distributions 
from transport models ( AMPT [TT], Hijing [12], UrQMD [13]) and a thermal model 
(Therminator [H]). By using several models with different physics implemented, we can 
study the effects of physics correlations and backgrounds that are trivially present in the 
data and that might modify purely statistical emission patterns, hke resonance decays, 
jet-production (Hijing), coalescence mechanism of particle production (AMPT), thermal 
particle production (Therminator), and hadronic rescatterring (AMPT, UrQMD). 

The kinetic coverage of protons and anti-protons used in our analysis is 0.4 < < 
0.8 GeV/c and \y\ < 0.5. In FiglH typical net-proton distributions for three centralities, 
0-5%, 30-40%, and 70-80%, of Au+Au collisions at ^/s]^ = 200 GeV are calculated 
with the UrQMD model. The shapes of net-proton distributions are different for the 
three centralities. For the most central collisions, — 5%, the net-proton distribution 
is wider compared to more peripheral collisions. The discrepancies in shapes will be 
reflected in the values of the different moments. 

Fig. 2 shows the < Npart > dependence of four moments (M, a, S, k) extracted 
from net-proton distributions of Au-|-Au collisions at y/s^N = 200 GeV for the various 
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Figure 1: Typical event by event net-proton multiplicity distributions of various 
centralities for Au+Au collisions at a/s/vat = 200 GeV calculated by UrQMD model. 
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Figure 2: Centrality dependence of various 
moments of Ap distributions for Au+Au 
collisions at -s/snn = 200 GeV from various 
models. The dashed lines represent the 
expectations for statistical emission. 
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Figure 3: Centrality dependence of various 
moments of Ap distributions for Au+Au 
collisions at various energies from the 
AMPT model. The dashed lines represent 
the expectations for statistical emission. 



models. M and a show a monotonic increase with < Npart > for all of the models, while 
iS" and K, decrease monotonically. In Fig. 3, we choose the default AMPT model to 
evaluate the centrality evolution of the various moments of net-proton distributions for 
various energies. M shows a linear increase with < Npart > and a decrease with ^/SNN■ 
a increases monotonically with < Npart > while it has non-monotonic dependence on 
^/sNN. S is positive and decreases with increasing < Npart > and ^/SNN■ The net- 
proton distributions become more symmetric for central collision and higher energies, k 
decreases with < Npart > and is similar for all energies. The dashed lines in Fig. 2 and 
Fig. 3 are resulting from Equs. (6)-(9) to evaluate the centrality evolution of the various 
moments. To apply our formulas, we fit the normalized mean value in Equ. (5) with the 
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function /(< A^^^^ >)• For AMPT String Melting (SM) [H], Hijing [12] and UrQMD 
[13] models, a 2nd order polynomial, /(< N^art >) = a < Npart >^ +b < N^art > is 
applied, while a linear function /(< Npart >) = a* < Npart > is employed for AMPT 
default [H] and Therminator [H] models. Once the function /(< Npart >) is obtained, 
the centrality evolution of the other moments is completely determined by Equs. (7)- 
(9). It is obvious that the centrality evolution of the various moments of net-proton 
distributions in Fig. 2 and in Fig. 3 can be well described by the dashed lines. 
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Figure 4: The kct^ of net-proton distribu- 
tions for Au+Au 7.7 GeV collisions ClS db 
function of < Npart > from various models. 
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Figure 5: The ko"^ of net-proton distribu- 
tions in Au+Au collisions function of 
^/SNN for various models. 



The KO"^ of the net-proton distributions of Au+Au collisions at ^/sNN = 7.7 GeV 
as a function of < Npart > is shown in Fig. 4 for various models, na"^ is constant 
with respect to < Npart > within the errors, which is consistent with the expectation 
from the IIES assumption. Fig. 5 shows the energy dependence of the na"^ of net-proton 
distributions for various models. The top of the figure shows the fiB value corresponding 
to the various energies. The values shown are averaged within the centrality range 
studied. The results from various models show no dependence on energy and are close 
to unity. This suggests that the kct^ of net-proton distributions is not affected very 
much by the non-QCP physics at different beam energies, such as the change of /^^ [1], 
and the collective expansion [15]. Note that the result from the pure thermal model, 
Therminator, is much closer to unity compared to others. Actually, if proton and 
anti-proton have independent poisson distributions, the difference of protons and anti- 
protons should distribute as a Skellam distribution [16] , for which na"^ is unity. A large 
deviation from constant as a function of < Npart > and collision energy for na"^ may 
indicate new physics, such as critical fluctuations. 
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4. Summary and Outlook 

Higher moments of the distribution of conserved quantities are predicted to be sensitive 
to the correlation length at QCP and to be related to the susceptibihties computed in 
Lattice QCD. Various non-QCP models (AMPT, Hijing, Therminator, UrQMD) have 
been applied to study the non-QCP physics background effects on the high moments 
of net-proton distributions. The centrality evolution of the high moments from models 
can be well described by the scaling derived from the IIES assumption and the moment 
products Sa, na/S and na^ of net-proton distributions are constant with respect to 
< Npart >■ /^o"^ is also found to be constant as a function of energy for various models. 

Our model study can serve as a background study of the behavior expected from 
known physics effects for the RHIC beam energy scan, that will span values of from 
100 to about 550 MeV. The presence of a critical point in that region may result in 
non-gaussian fluctuations and in correlated emission. Then the IIES assumption will 
break down. This is expected to lead to non-monotonic behavior of the observables 
studied here as a function of collision energy. 

Acknowledgments 

This work was supported in part by the U.S. Department of Energy under Contract 
No. DE-AC03-76SF00098 and National Natural Foundation of China under Grant No. 
10835005. BM thanks the Department of Atomic Energy, Government of India for 
financial support. 

References 

[1] Adams J et al, 2005 Nucl. Phys. A 757 102; Luo X F al, 2009 Physics Letters B 673 268; 

Mohanty B, 2009 Nucl. Phys. A 830 899C. 
[2] Aoki Y et at, 2006 Nature 443 675; Cheng M 2006 Phys. Rev. D 74 054507. 
[3] Ejiri S, 2008 Phys. Rev. D 78 074507; Bowman, E S, Kapusta J I, 2009 Phys. Rev. D 79 015202. 
[4] Stephanov M, 2005 Int. J. Mod. Phys A 20 4387; Fodor Z et al, 2004 JHEP 0404 50; Gavai R 

V, Gupta S, 2008 Phys. Rev. D 78 114503; de Forcrand Philippe and Philipsen O. 2002 Nucl. 

Phys. B 642 290. 

[5] Abelev B. I. et al., arXiv:0909.4131; STAR Internal Note - SN0493, 2009. 

[6] Stephanov M A ei al, 1999 Phys. Rev. D 60 114028; Berdnikov B and Rajagopal K, 2000 Phys. 
Rev. D 65 105017. 

[7] Stephanov M A, 2009 Phys. Rev. Lett. 102 032301. 

[8] Cheng M et al, 2009 Phys. Rev. D 79 074505. 

[9] Hatta Y and Stephanov M A, 2003 Phys. Rev. Lett. 91 102003. 
[10] Asakawa M, arXiv:0904.2089. 
[11] Lin ZW et al, 2005 Phys. Rev. C 72 064901. 

[12] Gyulassy M and Wang X N, 1994 Comput. Phys. Commun. 83 307. 

[13] Petersen H et al, arXiv:0805.0567. 

[14] Kisiel A et al, 2006 Comput. Phys. Commun. 174 669. 

[15] Abelev B. I. et al (STAR Collaboration), 2008 Phys. Rev. C 77 54901. 

[16] Skellam J G, 1946 Journal of the Royal Statistical Society 109 296. 



